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ABSTRACT 

M3 76 

Many of the satellites sent into orbit during the last few 
years are spinning for stabilization and or other reasons. 
The purpose of this paper is to determine the spin axis, a 
unit vector S°, of the satellite using ground tracking informa- 
tion. It will be shown that range rate information is adequate 
for this purpose. The determination of S° is necessary 
either to check the onboard sensors once the spacecraft is 
in orbit or in case of a slight onboard malfunction, the spin 
axis could not be determined otherwise. 
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DETERMINATION OF THE SPIN AXIS 
OF A SPINNING SATELLITE 


INTRODUCTION 

On many occasions it is most desirable or even necessary to determine the 
spin axis of a rotating satellite in space. The reasons for this are for instance 
check to onboard sensors, or to determine the spin axis in case of a malfunction. 
In this case many of the carried experiments can still be evaluated as long as 
the spin axis is known. 

One of the necessary conditions to do so with the method outlined here is 
that the radiating antenna element is displaced from the spin axis center since 
the Doppler modulation imposed and the range rate by the rotating satellite is 
utilized for the analysis. The range rate or Doppler as commonly used of a 
satellite in orbit can be determined with great precision. (See References 1, 2, 
3, and 4.) Errors in the order of cm/s can and have been achieved with cw- 
ranging systems. This fact will be utilized to determine the spin axis (unit 
vector S°) of a rotating satellite. 

The only assumptions made here are (a) that the antenna of the satellite is 
displaced from its axis of symmetry, (b) that the spacecraft is spinning only, 
and (c) that the satellite orbit can or has been determined. 


I. RANGE RATE MODULATION DUE TO SATELLITE SPIN 

Some of our satellites are spinning in space and have their radiating ele- 
ment, the antenna, displaced from its spin axis.* 

It will now be shown in this paper how the range rate r is modulated by the 
spinning motion of the spacecraft. Figure 1 and la depict the vector diagram of the 
antenna motion in detail as well as the angle \p = arc cos (r°, s° ) between the 
satellite spin axis s° and the satellite position vector 7 (of r°). This 
angle, as shown later, will be used for the determination of s°. Using vector 


*See reference 5 for centered turn stile antenna. 
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Antenna Motion Due to Spin 
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Figure la— Satellite and Antenna Details 
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notation , the antenna position vector d can be written as 


d = d Q cos ojAt + d(s°x d°) sin cu At 


( 1 ) 


where d° = d-d° * is the starting position of the spacecraft antenna at t = t # 
or At = 0. Equation 1 follows directly from Figure 1 as can be seen. Project- 
ing d onto the range vector r results in: 


Ar = (d- r°) =1 (d-r) 


( 2 ) 


Since, as assumed, the spacecraft is spinning, the value Ar changes with 
time (also due to orbital motion which is assumed to be negligible as is shown 
later), that is: 


=-L jr [(d-r) + (d-r)] - (d-r) fj 


(3) 


which is obtained by simply differentiating (2) with respect to time. 

Differentiating (1) also with respect to time and introducing it into (3) yields 
after some manipulation: 


Ar = - o)d(d° 


r°) sin oiAt + ojd(d° • (7° x s°)) cos ojAt + 


(7) (d°-7) - (4) r (d°-7°) 


(4) 


unit vector will always be designated with a zero superscript, that is x° = x / | x 
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Inspecting (4) shows that as long as 


|r I , r 

<< OJ and — < < CO 

r r 


( 5 ) 


the third and the fourth term of (4) can be neglected.* 

— 4 

Since d° is an arbitrary starting vector (unit vector) it will be chosen so 

that 


(d° -7° ) = 0 


( 6 ) 


or 


d °o 


1 


sin \p 


(r x s ) 


Introducing (6) into (4) neglecting the terms just mentioned yields: 


Ar = (o) d sin \p) cos a; At 


(7) 


with 


v r 

— < < u), — < < u 

r r 


as the "Range Rate Modulation" due to the spinning motion of a satellite. 


* Example: | 7 | = |v| (neglect earth rotation) = 8000 m/s 

n = 9 revolutions per min., co = 77n/30 = 1 sec 
r = 400 km (slant range) 
than: | r |/r i 2.10 -2 sec -1 << 1 sec' 1 
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If the range rate modulation Af. from the i th ground tracking station, is 
given as shown in Figure 2. The angle 0 . between the spin axis and the position 
vector r. (or r° ) can be determined from (7) that is 


0. = arc sin 



( 8 ) 


since a> = 277 /T is given by the same process. A better determination of co can 
be obtained from the power spectrum of the range rate information as shown in 
Figure 3. 

This angle 0. will now be used to determine the components of the unit 
vector s°, the satellite spin axis. 


II. DETERMINATION OF THE COMPONENTS OF THE SATELLITE 
SPIN AXIS UNIT VECTOR s° AND THEIR ERRORS 

The angle 0 . is certainly related to the satellite spin axis as can be seen 
from Figure 4. In reality, the mathematical relation can simply be expressed 
as a vector dot product, that is 


cos 0. = (r-° • s°) (9) 

where r° is the unit satellite position vector of the i th ground station. From 
(8), the angle 0 . can be determined, thus only r s° is unknown in equation (9). A 
straight forward solution of (9) would be to measure from 3 stations the Ar. and 
determine the components (Sj, s 2 , s 3 ) from (9). 

A much better solution can be obtained however by utilization of all informa- 
tion which is usually available, that is using much more than three measure- 
ments. In practice many of our ground stations take a large number (several 
hundreds) of range rate measurements. This in turn dictated a least square 
solution since more than the necessary three equations (measurements) are thus 
available. 

The dot product of (9) can also be written in component form as: 

cos = r n s i + r i2 s 2 + r i3 s 3 (10) 
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(b) Range Rate Modulation (Orbital Motion Subtracted) 
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Figure 2 
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Frequency for One Revolution 
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Figure 3-Power Density Spectrum 
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where of course the values r. ,, r r and s. , s. , s , are the components of 

l 1 1 2 *3 1 2 j 

the unit vectors r° and 1° respectively. 

Using as mentioned many measurements of Ar. and co { from which 4> i can 
be calculated using (8) one has a set of linear equations at hand. From (10) one 
obtains for n measurements of A f. and the following: 

a j = cos = r n s x + r 12 s 2 + r , 3 s 3 

a 2 = cos - r 21 Sj + r 22 s 2 + r 23 s 3 


a = cos 0 

n 'n 





(ii) 


Equation (11) suggested by itself a matrix treatment which will be followed 
from now on. In matrix form (11) reads simply 


(nx 1) 


R •«? 

( nX 3 ) (3x1) 


( 12 ) 


where the matrices A, R and S (the unit vector of the satellite spin axis) are as 
follows : 


a l 


r u 

r i2 

r 1 3 

a 2 


r 21 

r 22 

r 23 

* 

, R = 

* 

• 

* 

a 

L n _j 


_ r nl 

r n2 

r n3 


'(3x1) 


10 



In order to solve (12) in the "least square" sense, standard procedures are 
being applied. (See References 6, 7, 8, 9, 10, 11, and 12.) 

Since in practice (12) as such can and will never exist due to inhanced errors 
in the measurements as well as in the analyses one has to write instead: 

\nxl) ~ ®( n x3)’^(3xl) + £ (nxl) (13) 


where e is an error matrix with the following properties; the expectation 

E(e) = 0 (14) 


and the dispersion matrix 


V(s) 


(n Xn ) 


= E(€e T ) =a 2 . I 


(nXn ) 


(15) 


where I, v , is a unit matrix and a 2 , is the variance of the measurements 
made. This means in view of (11) to each equation in (11) an error term has to 
be added to encompass reality. 

Equation (13) is that equation which has to be used to determine the matrix 
S or the unit vector s° . The principle of "least squares" requires to solve for 
S in such a way that the sum of the errors squared (e T e ) is a minimum. ( e T is 
the transpose of c .) 

A necessary condition for this is that 

C' M - ° (16) 


or using (13) one obtains from (16) 


— [(A T - S t R t ) • (A - RS)] =0 
3S 


which yields after some calculations the "best” estimate for S , that is 


^(3x1) “ (® T ^(3 x 3) ^(3xn) ^(nxl) 


(17) 
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Equation (17) represents the solution to the problem of the determination of the 
satellites spin axis. 

In order to estimate the errors made in the determination of s°, its disper- 
sion matrix has to be calculated. From (17) and (15) one obtains: 

S = (R T R) _1 R t (RS + e) (18) 


Using (15) and introducing (18) one obtains after some calculations, 

V(S) = E(ee T ) = a 2 i (R T R) _1 < 19 ) 


for the dispersion matrix of S. 


HI. ESTIMATION OF THE MEASUREMENT VARIANCE 

In order to evaluate the dispersion matrix V(S) (3xl) the measurement 
variance has also to be estimated. 

From (9) it is to be remembered that Ar x and oj. are "measured" quantities 
which in turn have errors associated with them. These can be estimated from 
the noise of the data and this is to be assumed here as known. From equations 
(8) and (11) one can write 


sin i/i. 


1 

1 a i max 

d co. 


( 20 ) 


and 


cos ill. = a. 

r i 1 


Varying both equations, one obtains for the variation in a. 


8a. = 


Ar. 


/ Mr. V 

(dco,) 2 / 1 - [ — ax 

dco. 


8 ( Ar. ) — Ar. ( — - 

' i max' « max \ 


( 21 ) 


co. 

1 
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Assume now that the errors are uncorrelated and replace the value Ar. max by 
(Ar. + crAf. ) since it can never be zero due to the inherent errors in its 

' i max x max' 

determination, one obtains from (21): 


cr . 

ai 


Ar. + crAr, 

i max i max 



crAr. \ 2 

i max \ 


CO. 


crAr 2 + Ar 2 

i max i max 



1/2 


( 22 ) 


With the "Measurement" variance cr 2 . and the matrix R known the disposition 
matrix V(S) given by (19) can be evaluated. 

Equation (22) shows also the limitations for cr . ; (Ar. + crAf. / cod) has 
to be smaller than 1 to prevent too large errors. 


IV. CONCLUSIONS 

It has been shown how the range rate modulation of a spinning satellite with 
a displaced antenna can be used to determine its spin axis (unit vector) in space. 
An estimate of the spin axis errors associated with this method is also presented. 

It is interesting to note that no additional equipment on the satellite as well 
as on the ground is necessary. This method may well be suited for spin axis 
determination of communications and navigation satellites. It further could be 
used to check the axis of spin stabilized kick stages if deemed necessary by 
adding a Doppler transponder and antenna to the stage. 
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